An alternative calculation of the Casimir forces between birefringent plates 
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Barash has calculated the Casimir forces between parallel birefringent plates with optical axes parallel to the 
plate boundaries [Izv. Vyssh. Uchebn. Zaved., Radiofiz., 12, 1637 (1978)]. The interesting new feature of the 
solution compared to the case of isotropic plates is the existence of a Casimir torque which acts to line up the 
optical axes if they are not parallel or perpendicular. The forces were found from a calculation of the Helmholtz 
free energy of the electromagnetic field. Given the length of the calculations in this problem and hopes of an 
experimental measurement of the torque, it is important to check the results for the Casimir forces by a different 
method. We provide this check by calculating the electromagnetic stress tensor between the plates and showing 
that the resulting forces are in agreement with those found by Barash. 

PACS numbers: 12.20.Ds, 42.50.Lc, 46.55.+d 



I. INTRODUCTION 

Not long after Casimir's landmark paper (T) on the attrac- 
tive force between two parallel perfect mirrors, Lifshitz gener- 
alized the result to the case of isotropic dielectric plates J20. 
Lifshitz's formidable calculation was subsequently extended 
to allow for a third dielectric between the plates [4|, the case 
in mind being that of parallel plates immersed in a fluid. An 
obvious further generalization of the problem is to allow for 
anisotropic permittivities in the plates. The first exact solution 
for anisotropic plates was found in 1978 by Barash Q, who 
considered uniaxial (birefringent) plates with optical axes par- 
allel to the plate boundaries, separated by an isotropic dielec- 
tric. An expression was found for the Helmholtz free energy 
of the electromagnetic field at finite temperature, which was 
a function of the angle between the optical axes of the plates 
as well as of the plate separation. This means that there is a 
torque on the plates as well as a perpendicular force. These 
forces persist even at zero temperature due to the zero point 
energy of the quantum vacuum. The anisotropy thus leads to a 
new phenomenon in Casimir forces, a torque that acts to align 
the optical axes of the plates if they are not parallel or per- 
pendicular [5 6|. These results have been used to propose an 
experimental measurement of the Casimir torque using bire- 
fringent crystals immersed in ethanol (6J- 

The degree of complexity of this problem is considerable. 
This is reflected in the very lengthy expression found in 1 5 1 for 
the free energy, which we will not reproduce here (it can also 
be found in |6|, but note an important misprint lfl8l ). Given 
the length of the calculations and plans for an experiment to 
test the theory, it is highly desirable to have an independent 
solution of the problem using a different method. To date 
however, the only other analysis of this problem was a very 
simplified calculation [7| which gave a rough estimate of the 
Casimir forces. (The simpler problem of one isotropic and 
one birefringent plate with optical axis perpendicular to the 
boundary was analyzed recently in [8|.) 

In this paper we solve the problem exactly using Lifshitz's 
approach |2l[3j|4l, in which the electromagnetic stress tensor 
is found by a Green-function method. At room temperature, 



the forces for a realistic experimental set-up were found in O 
to be determined by the zero-temperature Casimir contribu- 
tion, with contributions from thermal radiation being negli- 
gible. We therefore ignore thermal effects and find the ex- 
act solution at zero temperature. Our approach has consider- 
able methodological interest, since we present a much simpler, 
more physical route to the stress tensor than is found in the 
standard Lifshitz theory ||2][3]|4). We developed this formal- 
ism in an analysis of the Casimir forces on moving plates |9l : 
here its general usefulness is demonstrated, since the solution 
obtained in this paper would have required considerably more 
labour if we had taken the usual approach [2 3 4]. Our treat- 
ment also generalizes the problem to allow for a magnetic 
response in the isotropic medium between the plates. Un- 
til recently an influence of magnetic permeability on Casimir 
forces was considered of only theoretical interest, but the 
development of metamaterials with engineered magnetic re- 
sponses raises new possibilities ll8l HTJl fTTl fl2l . 

From the electromagnetic stress tensor between the plates 
we obtain an expression for the electromagnetic energy; when 
there is no magnetic response our result agrees with the low- 
temperature limit of Barash's formula. This agreement is 
highly non-trivial: the formula for the energy in both solutions 
requires integrations and the integrands in each case are com- 
pletely different; numerical evaluation of the integrals shows 
that the two expressions give the same number. 

Section [II] sets out the problem and notation in detail. In 
Sec. |III] we show how the electromagnetic stress tensor can 
be found from the Green tensor of the vector potential. The 



Green tensor is calculated in Sec. IV by a novel method, and 



in Sec.[V]we present the stress tensor and electromagnetic en- 
ergy. 



II. NOTATION AND GEOMETRY 

We consider the arrangement depicted in Fig.[T] Two bire- 
fringent plates lie in the yz-plane, separated by a distance a. 
The optical axis of Plate 1 lies at an angle 8 to the z-axis 
while the optical axis of Plate 2 is directed along this axis. Be- 
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FIG. 1: Two plates birefringent plates lie in the yz-plane with con- 
stant separation a. The lines on the plates represent the directions of 
the optical axes, which also lie in the yz-plane. The axis of Plate 2 
lies in the z-direction while that of Plate 1 is rotated relative to this 
by 8. Between the plates is an isotropic medium of permittivity e and 
permeability fi. 



tween the plates there is an isotropic material with frequency- 
dependent permittivity e and permeability ji. 

The permittivity along the optical axis is denoted by ey and 
the permitivity in directions perpendicular to the optical axis 
by e ± ; these are also frequency dependent but to reduce the 
notational clutter we will not make this explicit. It follows 
that the dielectric tensors |[T3l fl4l of Plates 1 and 2 are, re- 
spectively, 
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A = 




cos 9 sin ( 
— sin 9 cos ( 
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(2) 



(3) 



((3) pro- 



is a rotation matrix about the x-axis. The angle 9 i 
duces the orientation of Plate 1 portrayed in Fig.[T] 

Luckily, the principle axes of a uniaxial crystal are fre- 
quency independent 1131 so the permittivity tensors are given 
by ([T])-([3]) at all frequencies. This property is crucial in 
making the problem tractable, although this was not explic- 
itly pointed out in 0. It does not hold for biaxial crys- 
tals (different permittivities in all three principle directions) — 
"dispersion of the axes" lfl3l . 



HI. LIFSHITZ THEORY 

The Casimir forces on the plates are determined by the 
expectation value of the electromagnetic stress tensor in the 
medium between the plates (4]. This expectation value is 
computed for the zero-temperature ground state of the elec- 
tromagnetic field, and the remarkable fact is that it is non- 
zero even though there are no electromagnetic fields present. 
The origin of this stress is the quantum zero-point energy of 
the electromagnetic field, which is modified by the materi- 
als lfT31l . From the classical Lagrangian of macroscopic elec- 
tromagnetism 



L = I d 3 x i(D • E - H • B) 



(4) 



one easily obtains the expression for the stress tensor by the 
usual methods of field theory [17|. The required quantity is 
the expectation value of the quantum version of this stress ten- 
sor: 



<r =e (D ®E 
1 



Mo 



(HOB) 



l(e (D-E) + Mo - i (H-B)). 



(5) 



The methodology of Lifshitz theory is to compute the expec- 
tation values in the stress |5]) using Green tensors of the vector 
potential A (with the scalar potential set to zero) 0. In quan- 
tum electrodynamics the basic computational quantity is the 
Feynman propagator or Green tensor G F , given by [3] 



(TA(r,t)®A(r7)) 



ihfiQ 
'~2tT 



G F (r,i;r',0- (6) 



This expression can be used to calculate field correlation func- 
tions which in turn determine the stress |5]l (see below). At 
finite temperature it can be shown |3] 01, through use of the 
fluctuation-dissipation theorem among other things, that the 
required correlation functions have a simple relation to the 
retarded Green tensor. In the zero-temperature case, how- 
ever, one can obtain this result much more simply using 
the fundamental starting point without having to invoke 
the fluctuation-dissipation theorem. We briefly describe the 
derivation. 

A Green tensor for the vector-potential wave equation sat- 
isfies 



G(r,i,r',fO 



duj G(r,r', oS) e 



-iu){t—t ) 



(7) 



V x 



or 



Vx -— e(«) G(r,r» =l<5(r-r'), (8) 



and depends on the material boundary conditions. From the 
monochromatic Maxwell equations leading to ([8]) one deduces 
the following property of the monochromatic Green tensor in 



G(r,r',-w*) = G*(r,r / ) w). 



(9) 



A complete set of boundary conditions for ([8]) requires not 
only the spatial boundary conditions given by the plates, but 
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also conditions in time. One thus has the usual choice |[T6l 
of retarded or advanced boundary conditions, or something 
more complicated. The Feynman Green tensor (|6|l can be con- 
structed from the retarded (G R ) and advanced (G A ) Green 
tensors, which are easier to calculate. In the limit r — > r', 
t — ► t', which is required to construct the stress tensor <(5j 
from the matrix element ([5J, the relationship is ifTTl 



limG F (r,i,r',i') 

t— *v 

= ilim(G ii (r,i,r',t') + G' 4 (r,t,r',0) 



(10) 



Since the boundary conditions are time symmetric, the ad- 
vanced Green tensor is found by time reversing the retarded 

case: 

G A (r,ty,t') = G R (v,~t,v',~t , ) (11) 
It follows from {TO]) and {FT} that 

limG F (r,i,r',i') = lim G K (r,t,r',t')> ( 12 ) 

so that only the retarded Green tensor is required to com- 
pute the stress. Once the retarded Green tensor is at hand 
the quadratic expectation values of the electric and magnetic 
fields in the stress tensor |5]) are found from [3 1 

(D(r)®E(r')) = -— f ^(i^G^r, r', (13) 
(H(r) ®B(r')) 

I" JO M«) 

These follow from (|6j, ( |T2| > and with a final Wick rotation 
to imaginary frequencies (uj = i£) in which |9]l is used. The 
switch to imaginary frequncies ensures that the integrals are 
well behaved El HH. 

IV. GREEN TENSOR 

The problem has now been reduced to that of finding the 
classical retarded Green tensor for the vector potential in the 
medium between the plates. With retarded boundary condi- 
tions, the solution of the monochromatic equation ([8]l has a 
simple physical meaning: an oscillating dipole at the point r' 
emits electromagnetic waves of frequency u and G R (r, r',ui) 
is the resulting vector potential at the point r. The second in- 
dex in G R represents the orientation of the dipole at r', while 
the first index represents the components of the vector poten- 
tial at r. Using this physical consideration it is clear from 
Fig.[T|that the solution will be a linear superposition of waves 
that have reflected off the plates, with the number of reflec- 
tions ranging from zero to infinity. To write down the solution 



we Fourier transform the Green tensor 

G R (x, x', u, V, i£) 

/>oo 

dy I dzG R {Y 1 r\iC) e- lu{ - v - yl) - lv(z - z,) 

(15) 



oo J — oo 



so that we decompose the waves emitted by the dipole into 
plane waves. In the absence of the plates the solution is the 
bare Green tensor 13] 01 



Gf (x,x ,u,v,i£) 



x>x > 
(x-x')g_^ x<x i 



G± 



1 
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(16) 



(17) 



(18) 



The two possibilities in (16i are linearly-polarized plane 
waves propagating to the right (first line) or to the left (sec- 
ond line), with wave vectors 



(±iw, u, v). 



(19) 



The imaginary x-component of the wave vectors is a conse- 
quence of the imaginary frequency, and in ( p~8] > we simply 
have the relation JtJiLu — ck. In physical terms the vac- 



uum solution ( 16 1 is trivial: it is the only way the dipole can 
propagate plane waves from x' to x. In the presence of the 



plates both plane waves in ( 16 1 will reflect off the plates and 
reverse direction, so the left-moving plane wave can propa- 
gate from x' to x even if x > x', with similar considerations 
applying to the right-moving plane wave. Consequently, both 
the right- and left-moving waves will appear in the solution 
regardless of whether x is greater or less than x', in contrast 



to the vacuum solution ( 16 1. Let R2 be the reflection operator 



(matrix) that transforms a right-moving plane wave at plate 2 
into the resulting reflected left-moving plane wave, and let R x 
be the reflection operator that transforms a left-moving plane 
wave at plate 1 into the reflected right-moving plane wave. We 
can now write down the solution (functional dependences are 
suppressed): 



G R =G R 



-w(x-x') j 



w(x-x') t 



+ (1 - e- 2a R!R 2 ) 1 

x (e- w[x - xl) g + + e- w{x+x ">n 1 g 

+ (l-e-^RaR!)^ 1 

x (e w( - x - x,) g- + e w( - x+x '~ 2a) R 2 G+ 



(20) 



The first line in ( 20 1 subtracts the left- or right-moving plane 
wave in ( 16 1, depending on whether x is greater or less than 



x . This subtraction is necessary because the direct propaga- 
tion, without reflections, of both the right- and left-moving 
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plane waves from x' to x is contained in the remaining terms 
in ( "20"} ; but only one of these propagations is possible, depend- 
ing on whether x is greater or less than x', and the first line 



The reflected wave at Plate2 has wave vector k_ (see ( 19 1), 



in (20 1 automatically subtracts the irrelevant one. The inverse 



matrices in ( 20 1 are geometrical series representing every pos- 
sible number of double reflections off both plates, the expo- 
nentials providing the propagation distance 2a for each double 
reflection. The initial right- and left-moving plane waves that 
leave x' reach x after both an even and odd number of reflec- 
tions; this explains the terms multiplying the inverse matrices 



in ( 20 1. Each term in ( 20 1, after the series expansion of the in- 



verse matrices, has an overall exponential factor that accounts 
for the propagation distance involved, with e~ ws , s > 0, rep- 
resenting a propagation distance s to the right for the initial 
right-moving plane wave, but to the left for the initial left- 
moving plane wave. 

The solution for the Green tensor, and thereby the Casimir 
forces, is thus determined once we compute the reflection op- 
erators Ri and R2, which describe the reflection of a linearly- 
polarized plane wave at each plate. This is a tedious business, 
but note that we have avoided the vastly more complicated 
task of having to solve the differential equation for the Green 
tensor with boundary conditions given by the plates. The use 



of physical reasoning to write down the solution ( 20 1 for the 
Green tensor is a novel feature of our approach to the Lifshitz 
theory for parallel plates. 



A. Reflection operator for Plate 2 

Consider first the reflection operator for Plate 2. There 
are two kinds of plane wave that can propagate in a bire- 
fringent crystal: the ordinary wave and the extraordinary 
wave (13\ [14*1 . When the plane wave with wave vector k + 
(see ( "H"") ) impinges on Plate 2 it produces an ordinary and an 
extraordinary wave in the plate with wave vectors lfl"3~l "141 



k 2o = (itfljo ,u,v), w 2o = VV + v 2 + e 2± K 2 , (21) 
k 2e = (iw 2e , u, v), w 2e = \Ju 2 + v 2 e 2ll /e 2± + e^K 
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(22) 



Although we are dealing with vector-potential waves, the 
electric field is proportional to A, so that A oscillates in the 
polarization direction. The vector potential of the ordinary 
and extraordinary waves can be written 



k-2o 



.4 



2o 



A 2e = A 



2<: 



-iw 2o I e 


iw 2e v 
uv 

v 2 + e 2 ±K 2 



ik2 -r+c^i 



'ik2 e -r+CKi 



(23) 



(24) 



which reveals their polarizations. Note that the ordinary wave 
is always polarized perpendicular to the optical axis. 



but its polarization is rotated relative to the incident wave. 
This rotation of the polarization by the plates is the reason 
why the final answer for the Casimir energy will be so much 
more complicated than in the isotropic case. 

To find the reflection operator for the k + plane wave we de- 
compose it into its component with polarization perpendicular 
to the plane of incidence (i?-polarization) and its component 
in this plane (S-polarization). The E- and S-polarization di- 
rections are given by the unit vectors 







n B2 = 



\Ju 2 + v 2 



1 



K^/eil ^Ju 2 + v 2 




(25) 



(26) 



The reason why there is no number subscript on is that 
this is also the £?-polarization direction for the k_ plane wave 
reflecting off Plate 1 . We must now solve the reflection prob- 
lem separately for the E- and £>-polarizations. This is done 
as in the isotropic case by imposing the standard boundary 
conditions that arise from the macroscopic Maxwell equa- 
tions ifTBI . From ("25") and ( |26] i we can take the incident E- 
and B-polarized waves to be 



A E = A E \ -v I e 

it 

,2 1 „,2\ 



ik_|_-r+c«t 



(27) 



i(u + v 2 ) 
vw 



Ar = An I uw I e lk +- r+CKt . (28) 



The answer for the reflection cannot be expressed in terms of 
a scalar reflection coefficient because of the rotation of the 
polarization referred to above. Nevertheless, for both inci- 
dent polarizations ( [27] ) and ( "28"" ), everything can be written in 
terms of two scalars which refer to the transmitted ordinary 
and extraordinary waves. The first of these scalars is the ratio 
A 2e /A 2o ; for the ^-polarization we denote this ratio by a 2 , 
and for the -polarization we denote it by j 2 : 



iu(ew 2o + e 2± w) 

a 2 = —, — 2 — . ^> 

v(ew^ + e 2 _ L w 2e w) 

ivw 2o (w + nw 2o ) 

e 2± un 2 (w + ^w 2e ) ' 



72 



(29) 
(30) 



The second scalar relates the transmitted ordinary wave to the 
incident wave; for the i?-polarization it is A 2o /A E , which we 
denote by (3 2 , and for the B-polarization it is ^20/^5, which 
we denote it by 5 2 : 
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02 = 



2ievw 

ew 2o w + e 2± u 2 + iuva 2 (ew + e 2± w 2e ) + e^i(v 2 + e 2± n 2 ) 

I 



52 = 



2ew(u 2 + v 2 ) 
euw 2o + iejivw^o + e 2± w(u + ij 2 vw 2e ) 



(3D 



The reflected waves for the two incident polarizations ( 27 1 and 
(|28]> are, respectively, 



-Vu 2 + v*(n E + m E2 )A E e *--+^, 
qi y/u 2 + v 2 (n B2 + m B2 )A B e &--H-c»rt 



where the vectors m E 2 and m B2 are 

j / e 2± /3 2 (u + ia 2 vw 2e )/e 
va.E2 = , „ = (3 2 (-iw 2o + a 2 uv) 
Vu 2 + v 2 y a 2 (3 2 (v 2 + c 2± k 2 ) 

j / £2±S 2 (-iu + j 2 vw 2e ) I 'e 

m B2 = — , , „ = -^2(^20 + «72™) 

Vu 2 + v 2 \ _ il2 s 2 ( v 2 + e2±K 2 ) 



(32) 
(33) 



(34) 



K-Jefi 



(35) 



The reflection operator R2 at Plate 2 for a general plane wave 
projects the wave to its components in the directions (|25] > and 
( po*} , reflects these components according to (|2"7}-(|35|>" and 
adds the reflected components; it is therefore given by 

R 2 = -(n E + m E2 ) ® n E - (n B2 + m B2 ) ® n B2 - (36) 



B. Reflection operator for Plate 1 

The reflection operator Tlx for Plate 1 can of course be ob- 
tained from R 2 . At Plate 1 the incident wave has wave vector 
k_ rather than k + , but this just means w — > —w compared to 
Plate 2. The i?-polarization direction of the incident wave is 



therefore still given by ( 25 1, but the i?-polarization direction 



risi is d26k with w — > — w 




(37) 



The only non-trivial issue is that the optical axis does not now 
lie in the z-direction. This is handled by changing the coordi- 
nate axes so that the optical axis lies along the new z-direction. 



The required passive coordinate transformation is a rotation 
about the x-axis by an angle 9; this transforms vector compo- 
nents with the matrix A -1 = A T , where A is given by IB). 
We denote quantities in the rotated basis by primes, so that the 
components of the incident wave vector k_ in this basis are 



k' = {—iw, u , v ), 



u cos a — v sin ( 



v = v cos 9 + u sin ( 



(38) 
(39) 



It is clear from (|2TT>— (|22]> that in the new frame the (left- 
moving) ordinary and extraordinary waves in Plate 1 have 
wave vectors 

k'io = (-iwi ,u',v'), wio = \Ju' 2 + v' 2 + ei ± n 2 , (40) 

k'le = (-iWi e ,u',v'), Wie = \J u' 2 + V ,2 £i ll /e 1± + eiyK 2 . 

(41) 

The reflection operator in the rotated frame, R' l5 is now ob- 
tained from R 2 by obvious replacements. We then find the 
reflection operator in the frame of Fig.[T]from 

Ri = ARiA T . (42) 

The result is 

Ri = -(n E + m.Ei) <g> n E - (n B i + m B1 ) ® n B1 , (43) 
where 



m E1 = 



m B1 = 



1 



y/u 2 + V 2 
1 



ei ± f3\{u' - iaxv'w le )/e 
Pi(iwi + a 2 u'v') 



, (44) 



nJTfi\/u 2 + v 2 



-6i(-w lo + ijiu'v') 
-i7i5i(u' 2 + ei ± K 2 ) 



(45) 

and at, (3%, 71, #1 are the Plate 1 versions of (|2"9)-(|3"T) in the 
rotated frame: 



ai 



iu'(ewi + ex ± w) 



v'(ew 2 lo 
2iev'w 



_w le w) ' 



7i 



ewioW + ei ± u' — iu'v'aii(ew + eij^Wie) + efj,(v' + ei ± K 2 ) 



iv'wi (w + fitvio) 

€l ± u'K 2 (w + HWi e ) 

Si 



(46) 



2ew(u' 2 + v' 2 ) 



eu'wio — it^\v'w\ + ei ± w(u' — ijiv'wi e ) 



(47) 
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With the reflection operators determined, we now have the 
solution for the Green tensor from ( |20| . Equation ( |20| con- 
tains inverse matrices and we evaluate these before using the 
Green tensor to compute the stress. This results in a lengthy 
expression for the Green tensor which we will not reproduce 
here. It is clear however from the form of the reflection op- 



erators (36 1 and (43 1, together with the geometrical-series ex- 



pansions of the inverse matrices in ( 20 1, that the Green tensor 
is largely constructed from scalar products of the vectors n^, 
iibi, ns2, m^i, m.E2, nisi and m^. Specifically, the fol- 
lowing scalar products can be usefully employed to expand 
the Green tensor: 





= (n B 4 


-m B1 ) • 


n B , 


(48) 


A 2 £_E 


= (n B 4 


- rn B2 ) • 


n B , 


(49) 




= (n B i 


+ m B i) 


• n B2 , 


(50) 


A 2 BB 


= (n S 2 


+ m B2 ) 


• n B i, 


(51) 


Mbe 


= m B i 


n E , 




(52) 


^2BE 


= m B2 


n E , 




(53) 



and these will appear in our solution for the Casimir stress and 
energy. 

V. STRESS TENSOR AND ENERGY 

Using the Green tensor we calculate the Casimir stress ten- 
sor er from Q an d (p"3j )-( 14 1. It is first necessary to drop the 
bare Green tensor (|16|> in ( 20 1, as this gives the diverging zero- 
point stress in the absence of the plates [3 |. The component 
a xx is the perpendicular Casimir force per unit area F on the 
plates; we find 



F 



he 

47T 3 



(Ik 



du dvw 



Ae 



law 



2B 



e 4aw _ Ae 2aw + B 



(54) 



where 



A — XieeMeE + AibbA 2 bb — 2XiBE^2BE, (55) 
B = (XiEE^iBB + X\ BE ) (A 2 _E£)A 2 _BB + A 2B£ ) . (56) 



If there is vacuum between the plates (e = fi = 1) the perpen- 
dicular force is always attractive. With a separating medium, 
however, even if [i = 1, the force can be attractive or repul- 
sive (6), as in the case of isotropic plates [4|. Note that the 
dependence of F on the plate separation a is entirely visible 



in ( 54 1, but the dependence on the orientation angle 9 is hid- 
den in the nested definitions that specify A and B. 



To compute the Casimir torque on the plates we require the 
Casimir energy E ; this is related to the pressure F and torque 
Qby 



F = 
Q = 



dE 

da ' 
dE 



(57) 
(58) 



We can therefore obtain E by integrating ( 54 1 with respect to 
a. This integration does not produce a function of integration 
depending on 9 but not on a, nor a constant of integration, 
since E must go to zero as a — > oo. Our final result is there- 
fore 



He f°° f°° f°° 
E =I^ dK du dv 

47r Jo J-oo J-oo 



2aw - - In (e iaw - Ae 2aw + B) 



(59) 



which determines the torque through (58 I. 

The previous exact analysis of this problem by Barash Q 
found a very different looking formula for the Casimir en- 
ergy. When /i = 1, the result (59 1 must agree with the zero- 
temperature limit of Barash's expression. This agreement can- 
not be seen analytically but we have verified it numerically. 



VI. CONCLUSIONS 

We have calculated the Casimir energy for two parallel bire- 
fringent plates with optical axes parallel to the boundary, al- 
lowing for the presence of an isotropic medium between the 
plates. This is only the second exact treatment of the problem, 
and our result is in agreement with the previous solution Q. 
Given the very different methods employed in the two cal- 
culations, this confirms the prediction of Casimir torque, an 
intriguing phenomenon that may be within reach of experi- 
ment J6). 
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